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in final states containing either two charged leptons and two jets, or one charged lep¬ 
ton, one neutrino and two jets, using proton-proton collision data at v^s = 8 TeV. 
The data, corresponding to an integrated luminosity of 19.7 fb^^, were recorded with 
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generation scalar leptoquarks with masses less than 1080 (760) GeV are excluded for 
/3 = 1.0 (0.5). Mass limits are also set for vector leptoquark production scenarios 
with anomalous vector couplings, and for R-parity violating supersymmetric scenar¬ 
ios of top squark pair production resulting in similar final-state signatures. These are 
the most stringent limits placed on the masses of vector leptoquarks and RPV top 
squarks to date. 
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1 Introduction 

The structure of the standard model (SM) of particle physics exhibits a symmetry between 
quarks and leptons. This paper reports on a search with the CMS detector at the CERN LHC for 
leptoquark (LQ) particles. These particles, which manifest a fundamental connection between 
quarks and leptons, are hypothesized by a variety of extensions to the SM such as grand uni¬ 
fied theories extended technicolor models HHHI, superstring-inspired models HJ], and 
composite models with lepton and quark substructure |T3|. Leptoquarks carry both baryon (B) 
and lepton (L) quantum numbers and thus couple to leptons and quarks. They carry fractional 
electric charge, are color triplets under SU(3)c, and can be either scalar or vector particles. 
Other properties such as their weak isospin, the helicity of the quarks and leptons to which 
they couple, and their fermion number F = (3B -|- L) depend on the specific structure of each 
model. Interpretations of direct searches for LQs at particle colliders rely on effective theories, 
such as the one described in Ref. HU, which require LQs to have renormalizable interactions, 
to obey SM gauge group symmetries, and to couple only to SM fermions and gauge bosons. In 
order to ensure proton stability, in effective theories LQs are generally constrained to conserve 
lepton and baryon numbers separately. Moreover, existing experimental limits IfT^ll^ on lep¬ 
ton number violation, flavor changing neutral currents, and other rare processes favor three 
generations of LQs with no intergenerational mixing, which is the scenario considered here. 

A search for pair production of first and second generation leptoquarks is performed in final 
states containing either two charged leptons and two jets, or one charged lepton, one neutrino 
and two jets, using proton-proton collision data at = 8 TeV. The data, corresponding to 
an integrated luminosity of 19.7 fb^^, were recorded with the CMS detector at the LHC. At 
hadron colliders, LQs would be produced in pairs or singly; this paper concentrates on LQ pair 
production. Recent CMS results for single LQ production are documented in Ref. IT/l . 

The production and decay of a scalar LQ are characterized by its mass (Mlq), its decay branch¬ 
ing fraction j6 into a charged lepton and a quark, and the Yukawa coupling characterizing 
the LQ-lepton-quark vertex. The interaction of scalar LQs with SM bosons is completely de¬ 
termined by these three parameters HU. The interaction of vector LQs with the SM bosons 
additionally depends on two anomalous couplings Ag and kq, which relate to the anomalous 
magnetic and electric quadrupole moments of the LQ that can be present in the gLQLQ and 
ggLQLQ vertices H8l , where g represents a gluon and LQ represents the anti-LQ. Four scenar¬ 
ios for the values of the anomalous couplings are typically considered in this case: minimal 
couplings (MC), A^ = 0, kg = 1; Yang-Mills (YM) type couplings, Ag = xg = 0; minimal- 
minimal (MM) couplings, Ac = xg = — 1; and the case of absolute minimal (AM) cross section 
with respect to the Ac, kq parameters for each value of LQ mass. 

LQ pair production arises predominantly through gluon-gluon fusion and quark-antiquark 
annihilation, shown in Fig. which have been calculated using next-to-leading order (NLO) 
QCD corrections |I19| . The dominant pair production mechanisms for scalar LQs do not depend 
on Af(j and the search sensitivity can be considered A^^-independent as long as A^^ is sufficiently 
large so that LQs decay within a few mm of the primary vertex. 

Other scenarios of physics beyond the SM could also lead to the prediction of particles with 
LQ-type couplings. One such theory is supersymmetry (SUSY), which postulates a symmetry 
between fermions and bosons, and predicts in some models the existence of quark superpart¬ 
ners (squarks), such as the top quark superpartner (top squark, t), decaying into LQ-like final 
states if R-parity is violated (RPV) l!20l . We consider one such model 11211 , where top squark 
decay is mediated by a Higgsino (H) with a mass Mg = Mj — 100 GeV with a 100% branching 
fraction. The Higgsino in turn produces an off-shell top squark, which decays to a charged 
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Figure 1: Dominant leading order diagrams for the pair production of scalar leptoquarks. 


lepton and a quark, as shown in Fig. The top squark decays via the RPV A-yj;, vertex, where 
A-yjj. represents the Yukawa coupling of the RPV term of the superpotential, and the ijk indices 
represent the family numbers of the interaction superfields, which correspond to l^e 

electron final state and A 232 for the muon final state. Limits have not previously been set on 
this model. 


t 



Figure 2: Diagram of the Higgsino-mediated top squark decay via the RPV Aj 32 (^=e) or A 232 
i£=fi) coupling. 


The final-state event signatures of the decay of pair-produced LQs can be classified as: dilep¬ 
ton and jets (both LQ and LQ decay into a charged lepton and a quark); single lepton, missing 
transverse momentum and jets (one LQ decays into a charged lepton and a quark, while the 
other decays into a neutrino and a quark); and missing transverse momentum and jets (both 
LQ and LQ decay into neutrinos and quarks). The three signatures correspond to branching 
fractions of 2j6(l — j 6 ), and (1 — j 6 )^, respectively. The charged leptons can be either elec¬ 
trons, muons, or tau leptons, corresponding to the three generations of LQs. Only final states 
containing electrons and muons are considered here, and two distinct signatures: one with two 
high transverse momentum (px) charged leptons and two high pj jets (denoted as ££jj), and 
the other with one high pj charged lepton, large missing transverse momentum, and two high 
Pt jets (denoted as ^vjj). These final states are analyzed in the context of scalar LQs, vector 
LQs 1221 and the RPV SUSY scenario previously mentioned. 

The most stringent limits on the pair-production of scalar LQs come from the LHC experi¬ 
ments. The ATLAS experiment excluded first (second) generation LQs with masses below 1050 
(1000) GeV for /3 = 1, and 900 (850) GeV for /3 = 0.5, using 20fb~^ of -y/s = 8 TeV data Il23| . 
Using = 7TeV proton-proton collisions data corresponding to an integrated luminosity 
of 5.0 fb^^, the GMS experiment excluded first- (second-)generation pair-produced scalar LQs 
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with masses below 830 (840) GeV for j6 = 1 and 640 (650) GeV for = 0.5 Il24l . GMS ex¬ 
cluded third-generation pair-produced scalar LQs with masses below 740 GeV for j6 = 1, using 
19.7 fb^^ of = 8TeV data Il25l . ATLAS excluded third-generation pair-produced scalar LQs 
with masses below 534 GeV for /I = 1, using 4.7 fb^^ of \/s = 7 TeV data 1261 . The HERA exper¬ 
iments HI Ii27l and ZEUS [28] produced A-dependent results for LQ models, excluding scalar 
LQ masses up to roughly 500-650 (300) GeV for A = 1.0 (0.3). Searches for scalar LQs have also 
been performed by the Tevatron experiments DO Il29] - l3ll and GDE Il33 - l34ll . The most stringent 
limits on vector LQs have been reported by DO Il35] - (37l and GDE Il38ll . 


2 The CMS detector 

The central feature of the GMS apparatus is a superconducting solenoid of 6 m internal diame¬ 
ter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip 
tracker, a lead tungstate crystal electromagnetic calorimeter (EGAL), and a brass and scintilla¬ 
tor hadron calorimeter (HGAL), each composed of a barrel and two endcap sections. Muons 
are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the 
solenoid. Extensive forward calorimetry complements the coverage provided by the barrel 
and endcap detectors. A more detailed description of the GMS detector, together with a def¬ 
inition of the coordinate system used and the relevant kinematic variables, can be found in 

Ref. |39l. 

The inner tracking system of GMS consists of a silicon pixel and strip tracker, providing the 
required granularity and precision for the reconstruction of vertices of charged particles in the 
range of the azimuthal angle 0 < (p < 2 k and pseudorapidity \ri\ < 2.5. The crystal EGAL 
and the brass and scintillator sampling HGAL are used to measure the energies of photons, 
electrons, and hadrons within \tj\ < 3.0. The electron momentum is estimated by combining 
the energy measurement in the EGAL with the momentum measurement in the tracker. The 
momentum resolution for electrons with pj^45 GeV from Z —> ee decays ranges from 1.7% for 
nonshowering electrons in the barrel region to 4.5% for showering electrons in the endcaps BOl . 

The GMS detector is nearly hermetic, which allows for a measurement of missing transverse 
momentum. The three muon systems surrounding the solenoid cover the region \i]\ < 2.4 
and are composed of drift tubes in the barrel region (|7/| < 1.2), cathode strip chambers in 
the endcaps (0.9 < |?7| < 2 .4), and resistive-plate chambers in both the barrel region and the 
endcaps {\ri\ < 1.6). Events are recorded based on a trigger decision using information from the 
GMS detector subsystems. The first level (LI) of the GMS trigger system, composed of custom 
hardware processors, uses information from the calorimeters and muon detectors to select the 
most interesting events in a fixed time interval of less than 4 }is. The high-level trigger (HLT) 
processor reduces the event rate from 100 kHz at LI to roughly 400 Hz. 


3 Data and simulation samples 

The data used in this paper correspond to an integrated luminosity of (19.7 ± 0.5) fb^^. The 
integrated luminosity is measured as described in Ref. [141 ]. 

Eor the searches in the eejj and ei/jj channels, events are selected by triggers requiring at least 
one electron with pj > 30 GeV, at least one jet with pj > 100 GeV, and at least one additional 
jet with Pj > 25 GeV. Eor the determination of the hadronic multijet background in the eejj and 
the evjj channels, events are selected using single-photon triggers, which require at least one 
EGAL energy deposit. 
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4 Event reconstruction and selection 


Events in the fiji]] and the f/vjj channels are selected if they pass a single-muon trigger selection 
that requires a muon with pj > 40 GeV, |p | < 2.1. There are no isolation requirements. This se¬ 
lection is also used to provide a sample of epjj events for the determination of the tf background 
in both ii]] channels. 

Simulated signal and background samples are produced and fully reconstructed using a sim¬ 
ulation of the CMS detector based on GEANT4 [42| . These simulations include additional 
collisions in a single bunch crossing (pileup) with a distribution matched to the number of 
pileup events observed during the various data-taking periods. 

Signal samples for scalar EQ masses from 300 to 1200 GeV in 50 GeV steps were generated at 
the leading-order (EO) level with the PYTHIA event generator Il43l and CTEQ6E1 |Ml parton 
distribution function (PDE) set. These samples are used to study the acceptance, while NEO 
cross sections are used for comparison in the limit-setting procedure. With the exception of 
the REV SUSY sample described below, PYTHIA 6.422 with the Z2 tune Il45l was used. The 
search limits being -independent, these samples were generated with a coupling strength 
= 0.3. The vector EQ signal samples were generated with the CalcHEP version 3.4 event 
generator H46ll and CTEQ6E PDE set using the model with vector EQ implemented in Ref. Il22l . 
Vector EQ masses between 200 and 1800 GeV were generated in 100 GeV steps, for the four 
scenarios of the anomalous couplings and kq described in Section]^ Samples of RPV SUSY 

events were produced with eejj and final state signatures. These samples were produced 

for top squark masses from 300 to 1000 GeV in 50 GeV steps using PYTHIA 8.175 071 and their 
decays were simulated withd MadGraph 5.1.1 Il48l . Top squark production cross sections 
were calculated at the NEO + next-to-leading-logarithm (NEE) level using PROSPINO Il4^ and 
the NLL-FAST program fSOHsTI, using the GTEQ6M PDE set. 

The main sources of background for these searches are tf, single top quark, Z/ 7 *+jets, W+jets, 
diboson (ZZ/WZ/WW)+jets, and multijet production. Backgrounds in the charmels from 
multijet production and tf events are estimated from data control regions, while single top 
quark, Z/ 7 *+jets, W+jets, and diboson (ZZ/WZ/WW)+jets backgrounds are estimated using 
simulated events. In the iv]] channel, the tf background is also estimated using simulated 
events. The simulated samples of tf, Z/ 7 *+jets, and W+jets are generated with MadGraph; 
single top quark samples (s-, i-, and tW- channels) are generated with POWHEG version 1.0 ||52l - 
[55l : and samples of VV, where V represents either a W or Z boson, are generated with PYTHIA. 
The simulations with MadGraph and PYTHIA use the GTEQ6E1 PDE set. The simulations 
with POWHEG use the GTEQ6M PDE set. 

The W+jets and Z/ 7 *+jets samples are normalized to next-to-NEO (NNEO) inclusive cross 
sections calculated with FEWZ version 3.1 1561 . Single top quark and VV samples are normal¬ 
ized to NEO inclusive cross sections calculated with mcfm version 6.6 Il57l - l60l . Results from 
Refs. IIM1l62l are used to normalize the tf sample at the NNEO + next-to-NEE level. 


4 Event reconstruction and selection 

Electron candidates are created by matching an electromagnetic cluster in the EGAE in t] and 
(pto a reconstructed track in the inner tracking system. The EGAE cluster must have a shower 
shape and longitudinal profile consistent with that of an electromagnetic shower. The matching 
reconstructed track can lack a hit in at most one pixel layer, and must be within 0.02 (0.05) cm 
of the matched primary vertex in the barrel (endcap). The resulting electron candidates are 
required to pass a set of criteria optimized for electrons with energies of hundreds of GeV BOl . 
In particular, they must have transverse momenta pj > 35 GeV and \r]\ < 2.5, excluding the 
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transition region between the barrel and endcap detectors, 1.442 < \r]\ < 1.560, where the 
electron reconstruction is suboptimal. The transverse momentum sum of tracks in a cone of 
AR = -y/ (A(/))^ + {At])^ = 0.5 around the electron candidate's track must be less than 5GeV, 
which reduces the chance of jets being misidentified as electrons. Tracks used in this momen¬ 
tum sum, known as tracker isolation, must be within 0.2 cm of the z coordinate of the electron 
candidate's matching primary vertex to eliminate tracks coming from other proton-proton colli¬ 
sions in the same bunch crossing. The transverse energy sum of the calorimeter energy deposits 
falling in the AR = 0.5 cone is required to be less than about 3% of the candidate's transverse 
energy This energy sum, known as calorimeter isolation, has an extra contribution accounting 
for the average contribution of additional proton-proton collisions in the same bunch crossing. 

Muons are reconstructed as tracks combining hit segments in the muon system and hits in the 
inner tracking system Il63]l . Muons are required to have pj > 45 GeV and to be contained in the 
fiducial volume used for the HLT muon selection, \i]\ < 2.1. In addition, muons are required 
to satisfy a set of identification criteria optimized for high pj. They require at least one muon 
detector segment be included in the muon track fit, and segments in at least two muon stations 
be geometrically matched to a track in the inner tracking system. Isolated muons are selected 
by requiring that the sum of the transverse momenta of all tracks in the tracker in a cone of 
AR = 0.3 around the muon track (excluding the muon track itself), divided by the muon pj, is 
less than 0.1. To have a precise pj measurement and to suppress muons from decays in flight, at 
least 8 tracker layers with associated hits are required, and at least one hit in the pixel detector. 
To reject muons from cosmic rays, the transverse impact parameter with respect to the primary 
vertex is required to be less than 2 mm and the longitudinal distance of the tracker-only track 
with respect to the primary vertex is required to be less than 5 mm, where the primary vertex is 
defined as the reconstructed vertex for which the pj sum of the assigned tracks is largest Il 6 ^ . 

Events are reconstructed using a particle-flow (PF) algorithm Il65l l66l . which identifies and 
measures stable particles by combining information from all GMS sub-detectors. The missing 
transverse momentum vector is defined as the projection on the plane transverse to the 
beams of the negative vector of the momenta of all particles reconstructed with the PF algo¬ 
rithm in the event, and the missing transverse energy is defined as the magnitude of the 

^miss vector. Jets are reconstructed using the anti-kj ||67l algorithm with a distance param¬ 
eter of 0.5. The jet energy is calibrated using the pj balance of dijet and 7 -l-jet events in both 
data and simulation Il69l . The PF jet energy resolution is 15% at 10 GeV, 8 % at 100 GeV, and 
4% at 1 TeV, to be compared to about 40%, 12%, and 5% obtained when the calorimeters alone 
are used for jet clustering. The leading (sub-leading) jet is required to have pj > 125 (45) GeV. 
All jets are required to have | //1 < 2.4. Furthermore, only jets having a spatial separation from 
electron or muon candidates of AR > 0.3 are considered. 

4.1 The ££]] channel 

An initial selection is made to obtain events containing at least two charged lepton candidates 
(either two electrons or two muons) and at least two jets for this channel. The two highest pj 
leptons and the two highest pj jets are considered as the decay products from a pair of FQs. 
They must satisfy the identification criteria described above. Further, the invariant mass of the 
two leptons. Mu, is required to be larger than 50 GeV. Muons are required to be spatially sepa¬ 
rated from one another by AR > 0.3. The scalar sum of the transverse momenta of the selected 
final state leptons and jets in the event Sj = Pt(^i) + Pt(^2) + PtOi) + PT(j2) is required to be 
larger than 300 GeV. No charge requirement is placed on the leptons. After this initial selec¬ 
tion, the signal-to-background separation is optimized by maximizing S / \/S + B, where S and 
B represent numbers of signal and background events, respectively This is done by varying 
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4 Event reconstruction and selection 


cuts on certain kinematic variables, and selecting the combination of cuts with the maximum 
S/-\/S + B. Three variables are optimized for each LQ mass hypothesis in both channels: 
St; Mu, used to remove most of the contribution from the Z/ 7 *+jets background; and 
defined as the smaller of the two lepton-jet invariant masses, given the combination that mini¬ 
mizes the LQ — LQ invariant mass difference. 

The eejj and mi]] channels are optimized separately and the optimized thresholds are summa¬ 
rized in Tables[^and|^ respectively For the mass hypotheses beyond 1 TeV the same set of final 
selections as those for the 1 TeV mass hypothesis is used. 

Table 1: Optimized thresholds for different LQ mass hypotheses of the eejj signal. 



300 

350 

400 

450 

500 

550 

600 

LQ mass [GeV] 

650 700 750 

800 

850 

900 

950 

> 1000 

St [GeV] 

435 

485 

535 

595 

650 

715 

780 

850 

920 

1000 

1075 

1160 

1245 

1330 

1425 

Mee [GeV] 

110 

110 

115 

125 

130 

140 

145 

155 

160 

170 

175 

180 

190 

195 

205 

^m.n [QgV] 

50 

105 

160 

205 

250 

290 

325 

360 

390 

415 

435 

450 

465 

470 

475 


Table 2: Optimized thresholds for different LQ mass hypotheses of the Hji]] signal. 


LQ mass [GeV] 



300 

350 

400 

450 

500 

550 

600 

650 

700 

750 

800 

850 

900 

950 

> 1000 

St [GeV] 

380 

460 

540 

615 

685 

755 

820 

880 

935 

990 

1040 

1090 

1135 

1175 

1210 

M,,,, [GeV] 

100 

115 

125 

140 

150 

165 

175 

185 

195 

205 

215 

220 

230 

235 

245 

M™" [GeV] 

115 

115 

120 

135 

155 

180 

210 

250 

295 

345 

400 

465 

535 

610 

690 


4.2 The iv]] channel 

Events in this channel are selected to contain exactly one charged lepton (electron or muon), 
at least two jets, and > 55 GeV. Leptons and jets must meet the criteria described above. 
Events containing a second lepton (electron or muon) are vetoed for the iv]] selections. In ad¬ 
dition, in order to reject events with mis-reconstructed the angle in the transverse plane 
between the direction of the leading px jet and the pjE'®® vector, A(/)(p^^®®,ji) is required to be 
larger than 0.5. For the same reason, the electron or muon and the are required to be 
separated by A<^(p^^®®, £) > 0.8. In the evjj charmel, the angular separation AR between the 
electron and either of the jets is required to be larger than 0.7 in order to reduce the contamina¬ 
tion from QCD multijet background in that channel. Events are required to have Mj > 50 GeV, 
where Mj, the transverse mass of the charged lepton and undetected particles, is defined as 
Mt = v/2p|£“'®® (1 — cos A^), where pj is the lepton px and Acp is the difference in azimuthal 
angle between the charged lepton momentum direction and the vector. Lastly, events are 
selected to have Sx > 300 GeV, where the scalar transverse energy Sj is defined in this case to 
be Sx = px(^) + Ex'®® + PtO'i) + PtO'i)- 

After this initial selection, the following variables are used to optimize a final selection for 
each LQ mass hypothesis using the method described above: Mj; Sx; and M^, defined as the 
invariant mass of the lepton-jet pair that minimizes the difference in the Mx of the lepton-jet 
and E™‘®®-jet pairs. The evjj channel uses as an additional optimization variable. 

The optimized thresholds for the evjj and the pvjj channels are summarized in Tables]^ and 
respectively. Mass hypotheses beyond 950 (1000) GeV for the evjj (pvjj) channel use the same 
set of final selections as those for the 950 (1000) GeV mass hypothesis. 
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Table 3: Optimized thresholds for different LQ mass hypotheses of the evjj signal. 



300 

350 

400 

450 

500 

550 

LQ Mass [GeV] 

600 650 700 

750 

800 

850 

900 

>950 

St [GeV] 

495 

570 

645 

720 

800 

880 

960 

1040 

1120 

1205 

1290 

1375 

1460 

1545 

Mej [GeV] 

195 

250 

300 

355 

405 

455 

505 

555 

600 

645 

695 

740 

780 

825 

Mt [GeV] 

125 

150 

175 

200 

220 

240 

255 

270 

280 

290 

295 

300 

300 

300 

[GeV] 

90 

95 

100 

no 

115 

125 

135 

145 

155 

170 

180 

195 

210 

220 


Table 4: Optimized thresholds for different LQ mass hypotheses of the signal. 



300 

350 

400 

450 

500 

550 

600 

LQ Mass [GeV] 
650 700 

750 

800 

850 

900 

950 

>1000 

St [GeV] 

455 

540 

625 

715 

800 

890 

980 

1070 

1160 

1250 

1345 

1435 

1530 

1625 

1720 

[GeV] 

125 

150 

175 

200 

225 

250 

280 

305 

330 

355 

380 

410 

435 

465 

490 

Mt [GeV] 

155 

180 

205 

225 

245 

260 

275 

290 

300 

310 

315 

320 

320 

325 

320 


5 Background estimation 

The main SM processes that can mimic the LQ signal in the ££]] charmels are: processes that 
lead to the production of genuine dilepton events such as Z/ 7 *+jets, tf, and VV+jets; and pro¬ 
cesses which produce either 0 or 1 genuine leptons and at least one hadronic jet which leads 
to a mis-identified lepton such as multijet events, single t production, and W+jets. The con¬ 
tributions from single top quarks, VV+jets, and W+jets are estimated from simulation and are 
small once the full event selection is applied. The contribution from the principal background, 
Z/ 7 *+jets, is estimated with using simulated events normalized to the data in a control region, 
where the non-Z/ 7 *+jets backgrounds have been removed from the data control region using 
the identical selection in simulation. The Z/ 7 *+jets simulation is rescaled to agree with this 
modified data sample at the £££]] initial selection level within a Z boson enriched region of 70 
(80) < <110 (100) GeV for the electron (muon) charmel. The resulting correction factor 

is Rz = 0.97 ± 0.01 (stat) for eejj, and Rz = 0.92 ± 0.01 (stat) for /i/rjj. The contribution from 
tt events with two leptons of the same flavor is estimated from a data sample containing one 
electron and one muon. This data sample is dominated by tt processes, which are expected to 
yield e/ijj events with the same probability as (eejj + /i/rjj) events. The data sample is therefore 
reweighted to account for: the different branching fraction of the e/ijj final state, which is twice 
that of the eejj or ^/ijj final states; the differences in electron and muon identification and iso¬ 
lation efficiencies; and the differences in trigger efficiencies. This sample can then be used to 
estimate the contribution from the tt process in the £££]] channels for both the initial and final 
selections in all kinematic distributions. 

The multijet background in the eejj channel is determined from a data control region containing 
exactly two electron candidates that pass loosened identification criteria on the cluster shape 
and no isolation requirements, and at least two jets. Each electron candidate in this sample 
is weighted by the probability that an electron candidate passing such loosened requirements 
additionally passes all final electron requirements. This probability is measured as a function of 
Pt in three//regions (II < 1.442,1.56 < |//| < 2.00, and 2.00 <\rj\< 2.50), using a data sample 
dominated by multijet events, collected with a single-photon trigger and containing one and 
only one electron candidate and two or more jets. The contribution from multijet processes 
in the ppjj channel is determined using a multijet-enriched data sample of same-sign dimuon 
events with no muon isolation criteria imposed. The same-sign nonisolated data sample is 
reweighted according to a same-sign/opposite-sign ratio and an isolation acceptance factor 
calculated using simulation. After reweighting, the same-sign nonisolated data sample is used 
to predict the multijet contribution to the final ppjj selection, which is shown to be negligible. 

All final state distributions in the eejj and ppjj charmels of the background prediction and of 
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data, at the initial selection level, have been studied and show agreement within uncertain 


ties. The specific distributions of Sj and are shown in Fig. 
discussed in the next section, are not included in these plots. 


Systematic uncertainties. 
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19.7 (8 TeV) 




Figure 3: Distributions of Sj (left) and (right) at the initial selection level in the eejj (top) 
and (bottom) channels. "Other background" includes: diboson, W+jets, 7 +jets, and single 
top quark contributions. The horizontal lines on the data points show the variable bin width. 


The primary backgrounds that can mimic the LQ signal in the iv']] channels fall into three cat¬ 
egories: events with genuine W bosons such as those from W+jets, tf, single top quark pro¬ 
duction, and WW and WZ processes; events with misidentified leptons and misreconstructed 
in the final state caused mostly by the misidentification of jets as leptons in multijet pro¬ 
cesses; and events with Z bosons such as those from Z/ 7 *+jets and ZZ processes, where only 
one lepton passes the identification and selection requirements. 

The contributions from the leading backgrounds (W+jets and tf) are determined using simu¬ 
lated events normalized to the data in control regions. The signal-depleted region 70 < Mj < 
110 GeV is used to determine both the W+jets and the tf normalization factors using two mu¬ 
tually exclusive selections. Selecting events with fewer than four jets produces a sample en¬ 
hanced with W+jets, and selecting events with at least four jets produces a sample enhanced 
with tf events. The results of these two selections are used to derive normalization factors from 
the following set of equations: 



19.7fb ' (8 TeV) 


500 1000 1500 2000 2500 3000 

Sj" [GeV] 


CMS 


Data 

Z/y* + jets 
tt 

Other background 
Multijet 
LQ, M = 450 GeV, p = 1.0 
LQ, M = 650 GeV, p = 1.0 


Ni = + RwlVi.w + IVi.o 

1^2 = Rtt^2,tt + Rw1V2,W + ^ 2,0 

where N,, jj, N,,W/ and Nip are the number of events in data, W+jets, tf, and other back¬ 
grounds passing selection i. The solution of the system yields the following normalization 
factors for the fn/jj channel: R^l = 0.99 ± 0.02 (stat) and Rw = 0.95 ± 0.01 (stat). Similar 
factors are obtained for the evjj channel: Rjj = 0.97 ± 0.02 (stat) ± 0.01 (syst) and Rw = 
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0.85 ± 0.01 (stat) ± 0.01 (syst), where the systematic uncertainties are associated with the esti¬ 
mate of the multijet background in this particular channel. The value of Rw in the evjj channel 
is affected by the lower efficiency of the trigger used in selecting W+jets events. 

The multijet background in the ei/jj channel is determined from data, using the previously de¬ 
scribed probability that an electron candidate satisfying loosened requirements also passes the 
final electron requirements. The probability is used to weight a sample of events containing: 
exactly one electron candidate passing the loosened identification criteria, at least two jets, 
and large The contribution from multijet processes is determined in the jiV]] channel 

using a sample of muon-enriched multijet simulated events with no muon isolation condition 
imposed. In the multijet-enriched region with < 10 GeV, the muon-enriched multijet 

simulated events are reweighted to agree with data, and a muon isolation acceptance rate is 
calculated using the data as the number of events passing the isolation condition divided by 
the total number of events. After re weighting and an adjustment by the muon isolation accep¬ 
tance factor, the nonisolated muon-enriched multijet simulated events are used to estimate the 
multijet contribution passing the final selection, which is determined in the /n/jj channel to be 
negligible. 

The contribution from the remaining backgrounds (diboson, single top quark, and Z/ 7 *+jets) 
is small and is determined entirely from simulation. 

As with the eejj and ^f/jj channels, all final state distributions in the ei/jj and jiv]] charmels of 
the background prediction and of data, at the initial selection level, have been studied and 
also show agreement within uncertainties. The specific distributions of Sj and for these 
channels are shown in Fig. Systematic uncertainties, discussed in the next section, are not 
included in these plots. 


19.7fb ' (8 TeV) 
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tt 
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2000 
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Figure 4: Distributions of Sj (left), and M^j (right) at the initial selection level in the evjj (top) 
and fiV]] (bottom) channels. "Other background" includes: diboson, Z/ 7 *+jets, and single top 
quark contributions. The horizontal lines on the data points show the variable bin width. 
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6 Systematic uncertainties 


6 Systematic uncertainties 

The sources of systematic uncertainty considered in this analysis are described below. To deter¬ 
mine the final uncertainty in signal and background predictions, each quantity is individually 
varied within its uncertainty, and the entire analysis is repeated to find the change in the pre¬ 
dicted number of background and signal events. 

The dominant sources of systematic uncertainty are the jet energy scale and resolution un¬ 
certainties, which are estimated by assigning a pj- and /^-dependent uncertainty in jet energy 
corrections, as described in Ref. 153, and by varying the jet pj by the uncertainty. The uncer¬ 
tainties in jet energy resolution are assessed by modifying the px difference between the re¬ 
constructed and particle-level jets by an p-dependent value 163 ranging between 5% and 30% 
for most jets. Charged-lepton momentum scale and resolution uncertainties also introduce un¬ 
certainties in the overall event acceptance. A px'dependent muon momentum uncertainty of 
5%x(px)/ where pj is expressed in TeV, and a px'dependent muon momentum resolution un¬ 
certainty ranging between 1 and 4% are used Il63 . For electrons in the ECAL barrel and endcap 
region, an energy scale uncertainty of 2% IfTOl and an electron energy resolution uncertainty 
of 10% |40| are used. The effects of these uncertainties are assessed by modifying the electron 
momentum scale and resolution in the simulation according to these uncertainties. A 2% per 
muon uncertainty in the muon reconstruction, identification, and isolation requirements, and a 
1% per muon uncertainty in the muon HLT efficiency are assigned in the ppjj and pi/jj channels. 
An additional uncertainty is assigned for the ppjj and pi/jj channels because of the effect on the 
muon momentum determination of the uncertainty on the alignment of the muon system. In 
simulation, a (p modulation can be seen in the difference between the inverse of the muon 
momentum as determined by the tracker with that determined by the tracker plus the muon 
system. Corrections were derived, but produced minimal differences, so instead a small uncer¬ 
tainty is added to account for possible alignment effects. In the ivjj analyses, the uncertainty in 
the charged lepton and jet energy and momentum scales and resolutions are propagated to the 
measurement of 

Other important sources of systematic uncertainty are related to the modeling of the back¬ 
grounds in the simulation. The uncertainties in the Z/ 7 *+jets, W+jets, and tt background 
shapes are determined using simulated MadGraph samples for which the renormalization 
and factorization scales and matrix element to parton shower matching thresholds have been 
varied up and down by a factor of two. The uncertainty of the scale factors for the normal¬ 
ization of the Z/ 7 *+jets background is determined to be 1% in both ££jj channels. A similar 
uncertainty for the normalization of the W+jets background is determined to be 2%(1%) in the 
epjj (ppjj) channel. The scale factor for the normalization of the tf background is determined to 
have an uncertainty of 2% in the e+jj and pvjj channels. The scale factor for the normali z ation 
of the epjj sample used for the tt background estimate in the eejj channel is determined to have 
an 8% uncertainty. 

The estimate of the multijet background from data in the eejj (e+jj) channel has an uncertainty 
of 60% (30%). This uncertainty is assessed by probing the precision of the method used to 
measure this type of background on an independent data control sample. 

An uncertainty in the modeling of pileup is determined by re-weighting the MC events to 
match with a number of pileup events 6% larger or smaller than what is observed in data, and 
an uncertainty of 2.6% is assigned to the value of the integrated luminosity 1411 . 

Lastly, the uncertainty in the signal acceptance, background acceptance, and cross section due 
to the PDF choice is estimated for signal (background) to be: 2% (3%) in the eejj channel; 3% 
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(3-25%) in the ei/jj channel; 2% (2-12%) in the jijiH channel; and 2% (1-21%) in the f/i/jj channel, 
following the PDF4LHC procedure Il^l72ll . 

The systematic uncertainties for both signal and background are summarized in Table for all 
channels, corresponding to the final selection optimized for Mlq = 650 GeV, which is repre¬ 
sentative of other high mass LQ values. 

Table 5: Systematic uncertainties (in %) for signal (S) and background (B) in all channels for the 
A^lq = 650 GeV final selection. 


Systematic 

Uncertainties 

eejj 

S [%] B [%] 

S [%] B [%] 

eyjj 

S [%] B [%] 

m 

S [%] 

i 

B [%] 

Jet energy scale 

0.30 

0.52 

0.42 

0.14 

1.6 

2.2 

0.02 

1.9 

Electron energy scale 

0.97 

6.4 

— 

— 

2.8 

3.3 

— 

— 

Electron Reco/ID/Iso 

4.0 

<0.01 

— 

— 

2.0 

<0.01 

— 

— 

Muon momentum scale 

— 

— 

0.63 

1.7 

— 

— 

0.19 

13 

Muon Reco/lD/lso 

— 

— 

4 

0.48 

— 

— 

2.0 

0.19 

Jet resolution 

0.01 

0.23 

0.23 

0.86 

0.09 

0.46 

0.78 

2.2 

Electron resolution 

0.46 

0.22 

— 

— 

0.61 

0.53 

— 

— 

Muon resolution 

— 

— 

0.14 

0.39 

— 

— 

0.15 

7.1 

Muon alignment 

— 

— 

0.1 

0.54 

— 

— 

1.0 

2.8 

Trigger 

<0.01 

<0.01 

<0.01 

<0.01 

<0.01 

<0.01 

1.0 

0.10 

tt normalization 

— 

2.1 

— 

0.35 

— 

1.5 

— 

0.60 

tt shape 

— 

— 

— 

— 

— 

3.0 

— 

1.4 

W+jets normalization 

— 

<0.01 

— 

0.01 

— 

0.12 

— 

0.63 

W+jets shape 

— 

<0.01 

— 

0.23 

— 

0.87 

— 

13 

Z/y* +jets normalization 

— 

0.75 

— 

0.59 

— 

<0.01 

— 

0.07 

Z/ 7 *+jets shape 

— 

12 

— 

12 

— 

<0.01 

— 

1.5 

Multijet modeling 

— 

0.10 

— 

— 

— 

5 

— 

— 

PDF 

2.0 

2.1 

2.0 

4.8 

3.0 

13 

3.0 

5.1 

Pileup 

0.04 

0.38 

0.16 

0.22 

0.14 

1.2 

0.06 

1.3 

Integrated luminosity 

2.6 

0.10 

2.6 

0.31 

2.6 

0.47 

2.6 

0.25 

Total 

5.3 

14 

5.2 

13 

5.5 

15 

4.7 

21 


7 Results 

Data and background predictions are compared for every channel and each mass optimization 
point, after the optimized final selection criteria are applied to both signal and background. 
The first part of this section details such comparisons. There are no significant deviations from 
SM background predictions. Limits are set on the cross section times branching fraction for 
the hypothesis of scalar LQ pair production as a function of Mlq and fi. The expected and 
observed limits for scalar LQ pair production are detailed in the second part of this section. 
Additional interpretations of the results in the context of vector LQ pair production and of 
RPV SUSY production with ££]] and £i/jj signatures are described in the last part of this section. 

7.1 Data and background comparison 

Agreement is found between data and background predictions in both the and jiV]] chan¬ 
nels, as shown in Fig. j^for the fifi]] channel, which displays Sj and for signal LQ masses 
of 450 and 650 GeV, and in Fig. for the channel, which displays Sj and Mf,j for the same 
signal LQ mass points. 

The numbers of events selected in data, and the various backgrounds at final selection as a 
function of Mlq for the and f/vjj channels are summarized in Tables]^ andrespectively. 
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7 Results 






Figure 5: Distributions of St (left) and (right) for the final selection for a LQ mass of 

450 GeV (top) and 650 GeV (bottom) in the mi]] channel. The dark shaded region indicates the 
statistical and systematic uncertainty in the total background prediction. "Other background" 
includes diboson, W+jets, and single top quark contributions. The horizontal lines on the data 
points show the variable bin width. 
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19.7fb '(8TeV) 19.7 (8 TeV) 



Figure 6: Distributions of Sj (left) and (right) for the final selection for a LQ mass of 450 GeV 
(top) and 650 GeV (bottom) in the fiV]] channel. The dark shaded region indicates the statistical 
and systematic uncertainty in the total background prediction. "Other background" includes 
diboson, Z/ 7 *+jets, and single top quark contributions. The horizontal lines on the data points 
show the variable bin width. 
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7 Results 


Since mass hypotheses at 1 TeV and beyond share the same final selections, they also share the 
same background yields. 

Table 6: Event yields for the analysis for /3 = 1.0 for all values of Mlq. Uncertainties are 
Poisson uncertainties in the simulated background, except for the second uncertainty for "All 
background", which gives the total systematic uncertainty as detailed in Section]^ Systematic 
uncertainties are dominated by jet energy scale and simulation shape uncertainties. 


Mlq [GeV] 

Signal 

Z/ 7 *+jets 


tt 

VV, W, single t 

All background 

Data 

300 

16240 

± 

110 

819.0 

± 

9.2 

666 

± 19 

88 

± 5.4 

1573 

± 22 ± 56 

1659 

350 

7570 

± 

48 

351.7 

± 

6.0 

405 

± 15 

58.3 

± 4.5 

815 

± 17 ± 21 

797 

400 

3658 

± 

22 

200.4 

± 

4.5 

202 

± 11 

31.5 

± 3.3 

434 

± 12 ± 17 

439 

450 

1816 

± 

11 

110.2 

± 

3.3 

103.7 

± 27.5 

20.6 

± 2.7 

234.6 

± 8.6 ± 11.2 

233 

500 

938.1 

± 

5.5 

69.9 

± 

2.6 

61.0 

± 5.7 

13.2 

± 2.2 

144.2 

± 6.7 ± 8.5 

135 

550 

498.8 

± 

2.9 

39.6 

± 

1.9 

29.4 

± 3.9 

8.0 

± 1.8 

77 

± 4.7 ± 5.2 

84 

600 

274.7 

± 

1.6 

25.8 

± 

1.5 

14.8 

± 2.8 

6.5 

± 1.6 

47.1 

± 3.6 ± 4.5 

47 

650 

157.1 

± 

0.9 

17.1 

± 

1.2 

10.3 

± 2.3 

4.3 

+1.4 

-1.3 

31.7 

+2.9 \ A 'y 
-2.9 ^ 

25 

700 

89.49 

± 

0.52 

10.71 

± 

0.98 

7.0 

± 2.0 

2.9 

+1.2 

-1.0 

20.6 

tli ± 4.3 

15 

750 

52.39 

± 

0.30 

6.95 

± 

0.79 

2.20 

± 0.98 

1.1 

+0.8 

-0.56 

10.3 

+('4 ± 2.7 

11 

800 

31.3 

± 

0.18 

3.90 

± 

0.59 

1.08 

± 0.62 

0.77 

+0.73 

-0.46 

5.8 

+lJ ± 1-55 

9 

850 

18.99 

± 

0.11 

1.96 

± 

0.39 

0.0 

+0.65 

-0.0 

0.75 

+0.73 

-0.46 

2.71 

+0.97 4^ 1 Q7 
-0.65 ^ 

5 

900 

11.290 

± 

0.067 

1.10 

± 

0.29 

0.0 

+0.65 

-0.0 

0.30 

+0.60 

-0.20 

1.38 

-1% ± 0.44 

3 

950 

6.907 

± 

0.041 

0.76 

± 

0.25 

0.0 

+0.65 

-0.0 

0.12 

+0.78 

-0.12 

0.87 

+032 ± 0.49 

1 

1000 

4.175 

± 

0.026 

0.41 

± 

0.18 

0.0 

+0.65 

-0.0 

0.0 

+0.77 

-0.0 

0.41 

+0-91 + 0 27 
-0.18 ^ 

0 

1050 

2.778 

± 

0.017 

0.41 

± 

0.18 

0.0 

+0.65 

-0.0 

0.0 

+0.77 

-0.0 

0.41 

+0-91 + n 27 
-0.18 ^ 

0 

1100 

1.860 

± 

0.011 

0.41 

± 

0.18 

0.0 

+0.65 

-0.0 

0.0 

+0.77 

-0.0 

0.41 

+0-91 + n 27 
-0.18 ^ 

0 

1150 

1.2471 

± 

0.0072 

0.41 

± 

0.18 

0.0 

+0.65 

-0.0 

0.0 

+0.77 

-0.0 

0.41 

+0-91 + n 27 
-0.18 ^ 

0 

1200 

0.8202 

± 

0.0047 

0.41 

± 

0.18 

0.0 

+0.65 

-0.0 

0.0 

+0.77 

-0.0 

0.41 

+0-91 + n 27 
-0.18 ^ 

0 


In both the eejj and evjj channels, a broad data excess is observed for the selections optimized 
for a LQ mass greater than about 400 GeV, as shown in Figs. and for two chosen selec¬ 
tions, and in Tables]^ andThis excess is most significant in the selection optimized for a LQ 
mass of 650 GeV, where for the eejj (ei/jj) channel 20.5 ± 2.1 (stat) ± 2.8 (syst) (7.5 ±1.2 (stat) ± 
1.1 (syst)) events are expected and 36 (18) events are observed, with a significance of 2.3 (2.6) 
standard deviations. 

An investigation of the kinematic distributions in both channels shows that the excesses are 
background-like. In particular, unlike a LQ hypothesis, the excesses do not peak sharply in 
the and the Mq distributions, as shown in Fig. ^ For comparison, the distributions that 
would result from a LQ mass hypothesis of 650 Gev and j6 = 0.075 are also shown (this is 
the value of j6 that, for a LQ mass of 650 GeV, would produce 10 events in the evjj selection 
optimized for such a LQ mass, which is about the size of the excess). The intrinsic width of 

scalar LQs is ^ x Mlq. The LQ signal events were generated with = 0.3. This corresponds 
to an intrinsic width of about 1.2 GeV for a LQ with mass close to 650 GeV, which is negligible 
compared to the experimental resolution. Significantly higher values of (and consequently 
broader LQs) are strongly limited in this mass range by results from the HERA experiments Il27l 

Further investigations of the characteristics of the data that survives the selections optimized 
for a LQ mass of 650 GeV show that there are two events containing same-sign electrons out of 
the 36 events, and we expect the SM background to contribute about two events with same- 
sign electrons out of the about 20 predicted events, because of charge misidentification. We 







7.2 Exclusion limits on scalar LQ pair-production 


15 


Table 7: Event yields for the jiv]] analysis for /3 = 0.5 for all values of Mlq. Uncertainties are 
Poisson uncertainties in the simulated background, except for the second uncertainty for "All 
background", which gives the total systematic uncertainty as detailed in Section]^ Systematic 
uncertainties are dominated by the jet energy scale and simulation shape uncertainties. 


Mlq [GeV] 

Signal 

W+jets 

tt 

VV, Z, single t 

All background 

Data 

300 

5089 ± 

58 

1102 ± 

22 

1853 ± 15 

331.3 ± 8.3 

3286 ± 28 ± 185 

3549 

350 

2352 ± 

25 

472 ± 

14 

640.0 ± 8.5 

159.8 ± 5.7 

1272 ± 18 ± 70 

1451 

400 

1064 ± 

11 

213.9 ± 

9.6 

259.5 ± 5.4 

84.6 ± 4.3 

558 ± 12 ± 38 

668 

450 

526.7 ± 

5.5 

115.7 ± 

7.1 

116.3 ± 3.6 

44.8 ± 2.9 

276.7 ± 8.5 ± 22 

313 

500 

263.6 ± 

2.8 

66.4 ± 

5.3 

56.1 ± 2.5 

25.1 ± 2.1 

147.6 ± 6.2 ± 12.6 

173 

550 

142.7 ± 

1.5 

43.8 ± 

4.4 

26.1 ± 1.7 

14.3 ± 1.6 

84.3 ± 5.0 ± 9.4 

93 

600 

78.1 ± 

0.8 

20.3 ± 

2.7 

13.7 ± 1.2 

8.0 ± 1.1 

42 ± 3.2 ± 5.8 

57 

650 

44.62 ± 

0.46 

14.0 ± 

2.3 

7.97 ± 0.95 

4.34 ± 0.72 

26.3 ± 2.6 ± 5.2 

36 

700 

25.27 ± 

0.26 

9.1 ± 

1.8 

5.20 ± 0.76 

9 70 + 0.64 
- 0.46 

1 qq + 0.60 
_ o .40 

17 ± 2.0 ± 4.7 

25 

750 

15.04 ± 

0.15 

7.0 ± 

1.6 

2.82 ± 0.56 

11.7 ± 1.8 ± 5.1 

15 

800 

9.080 ± 

0.093 

4.5 ± 

1.4 

1.47 ± 0.41 

4-D4 _ o .37 

7.6 ± 1.5 ± 3.5 

11 

850 

5.493 ± 

0.056 

1.08 ± 

0.54 

1.04 ± 0.35 

1 +0.55 

4-40 _ Q 32 

3.28 ± 1.04 

7 

900 

3.370 ± 

0.035 

0.62 ± 

0.44 

0.92 ± 0.32 

f) q +0.53 
-0.29 

2.44 toll ± 0.89 

3 

950 

2.111 ± 

0.022 

0.4 ± 

0.4 

0.44 ± 0.22 

0 51 

- 0.21 

1.35 ± 0.6 

3 

1000 

1.322 ± 

0.014 

0.4 ± 

0.4 

0.26 ± 0.18 

0 51 

- 0.21 

1.17 +°|j ± 0.56 

3 

1050 

0.9338 ± 

0.0092 

0.4 ± 

0.4 

0.26 ± 0.18 

0 51 

- 0.21 

1.17 +°|j ± 0.56 

3 

1100 

0.6507 ± 

0.0062 

0.4 ± 

0.4 

0.26 ± 0.18 

0 51 

- 0.21 

1.17 +°|j ± 0.56 

3 

1150 

0.4457 ± 

0.0041 

0.4 ± 

0.4 

0.26 ± 0.18 

0 51 +0.49 
-0.21 

1.17 +°|j ± 0.56 

3 

1200 

0.3097 ± 

0.0028 

0.4 ± 

0.4 

0.26 ± 0.18 

n 51 

- 0.21 

1.17 tofi ± 0.56 

3 


have also verified that the excess is not enhanced if we require that the jets are identified as 
b-quark jets using the combined secondary vertex b-tagging algorithm Il7^ . 

A recently published search for heavy neutrinos and W bosons with right-handed couplings f7^ 
also observed an excess in the number of selected eejj events compared to the expectation 
from SM backgrounds. However, the excess in Ref. f/Sl is mostly localized in the region 
1.8 < Meejj < 2.2 TeV, where Mgejj is the invariant mass of the 2 leading electrons and 2 lead¬ 
ing jets, while the excess observed in this analysis with the selection optimized for LQ mass of 
650 GeV is broadly distributed between Mgejj values of 1 and 2 TeV. Furthermore, only 30% of 
the events populating the excess region in Ref. f74l survive the Mlq = 650 GeV selection. 

In summary, the kinematic properties of the data in the excess regions for the eejj and the 
ei/jj channels are not found to be consistent with a LQ signal, and the size of the data excess 
is significantly less than that expected for a LQ with a mass of 650 GeV and j6 > 0.5. In the 
following section, limits are set on LQ production for both first and second generation. 


7.2 Exclusion limits on scalar LQ pair-production 


Upper limits are set on the scalar LQ production cross sections a using the asymptotic GLs 
modified-frequentist approach |[^[^. A log-normal probability density function is used to 
integrate over the systematic uncertainties described in Section]^ Uncertainties of statistical na¬ 
ture are described with gamma distributions with widths determined by the number of events 
in signal and background simulated samples or observed in data control regions. 


The 95% confidence level (GL) upper limits on tr x or cr x 2j6(l — j}) as a function of LQ 
mass are shown together with the NLO predictions for the scalar LQ pair production cross 
section in Fig. 10 for the eejj and evjj channels, and in Fig. 11 for the fiji]] and /n/jj channels. 
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7 Results 


19.7fb '(8TeV) 19.7 fb“'(8 TeV) 






Figure 7: The Sj (left) and (right) distributions for events passing the eejj selection op¬ 
timized for Mlq = 450 GeV (top) and Mlq = 650 GeV (bottom). The dark shaded region 
indicates the statistical and systematic uncertainty in the background total prediction. "Other 
background" includes diboson, W+jets, and single top quark contributions. The horizontal 
lines on the data points show the variable bin width. 


Table 8: Event yields for the eejj analysis for j6 = 1.0 for all values of Mlq. Only statistical uncer¬ 
tainties are reported, except in the "All background" column, where systematic uncertainties 
are also reported. 


Mlq [GeV] 

Signal 

Z/ 7 *+jets 

tt 

Multijet 

VV, W, single t 

All background 

Data 

300 

13560 

± 

80 

462.2 

± 

7.4 

724 ± 20 

5.280 

± 

0.052 

62.1 

± 

4.6 

1254 ± 22 ± 76 

1244 

350 

6474 

± 

33 

332.1 

± 

6.2 

352 ± 14 

3.220 

± 

0.036 

37.7 

± 

3.6 

725 ± 16 ± 48 

736 

400 

3089 

± 

15 

203.2 

± 

4.8 

153.7 ± 9.1 

1.700 

± 

0.023 

23.8 

± 

2.9 

382 ± 11 ± 27 

389 

450 

1508 

± 

7.2 

112.9 

± 

3.5 

86.9 ± 6.9 

0.890 

± 

0.016 

11.8 

± 

2.0 

212 ± 8.0 ± 18 

233 

500 

767.4 

± 

3.6 

66.5 

± 

2.7 

47.2 ± 5.1 

0.490 

± 

0.011 

7.4 

± 

1.6 

122 ± 6.0 ± 9.3 

148 

550 

410.5 

± 

1.9 

37.4 

± 

2.1 

25.8 ± 3.7 

0.2800 

± 

0.0084 

3.7 

± 

1.1 

67.2 ± 4.4 ± 5.2 

81 

600 

225.7 

± 

1.0 

22.2 

± 

1.6 

14.2 ± 2.8 

0.1500 

± 

0.0065 

3.12 

± 

1.00 

39.7 ± 3.4 ± 3.3 

57 

650 

125.90 

± 

0.58 

14.0 

± 

1.2 

5.4 ± 1.7 

0.0760 

± 

0.0040 

1.05 

± 

0.47 

20.5 ± 2.1 ± 2.8 

36 

700 

72.88 

± 

0.33 

8.16 

± 

0.93 

4.3 ± 1.5 

0.0450 

± 

0.0029 

0.21 

± 

0.12 

12.7 ± 1.8 ± 2.3 

17 

750 

43.10 

± 

0.20 

4.88 

± 

0.69 

1.55 ± 0.90 

0.0260 

± 

0.0023 

0.078 

± 

0.038 

6.5 ± 1.1 ± 1.2 

12 

800 

26.17 

± 

0.12 

2.93 

± 

0.52 

1.04 ± 0.73 

0.0190 

± 

0.0022 

0.078 

± 

0.038 

4.06 ± 0.90 ± 0.93 

7 

850 

15.980 

± 

0.072 

2.34 

± 

0.48 

0.52 ± 0.52 

0.0110 

± 

0.0015 

0.042 

± 

0.028 

2.91 ± 0.71 ± 0.74 

5 

900 

9.813 

± 

0.044 

1.23 

± 

0.36 

0.52 ± 0.52 

0.0069 

± 

0.0012 

0.022 

± 

0.020 

1.77 ± 0.63 ± 0.39 

3 

950 

6.086 

± 

0.028 

0.89 

± 

0.29 

0.00 f 

0.00450 

± 

0.00085 

0.022 

± 

0.020 

0.91 ± 0.28 

1 

1000 

3.860 

± 

0.018 

0.56 

± 

0.22 

0.00 iji'* 

0.00370 

± 

0.00082 

0.0025 

± 

0.0025 

0.57 tlit ± 0.18 

1 

1050 

2.576 

± 

0.011 

0.56 

± 

0.22 

0.00 f J 

0.00370 

± 

0.00082 

0.0025 

± 

0.0025 

0.57 ± 0.18 

1 

1100 

1.6940 

± 

0.0072 

0.56 

± 

0.22 

0.00 

0.00370 

± 

0.00082 

0.0025 

± 

0.0025 

0.57 ± 0.18 

1 

1150 

1.1270 

± 

0.0047 

0.56 

± 

0.22 

0.00 f J i'* 

0.00370 

± 

0.00082 

0.0025 

± 

0.0025 

0.57 ± 0.18 

1 

1200 

0.7500 

± 

0.0030 

0.56 

± 

0.22 

0.00 f 

0.00370 

± 

0.00082 

0.0025 

± 

0.0025 

0.57 ± 0.18 

1 
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Figure 8: The Sj (left) and Mq (right) distributions for events passing the full evjj selection 
optimized for Mlq = 450 GeV (top) and Mlq = 650 GeV (bottom). The dark shaded region 
indicates the statistical and systematic uncertainty in the total background prediction. "Other 
background" includes diboson, Z/ 7 *-i-jets, and single top quark contributions. The horizontal 
lines on the data points show the variable bin width. 


Table 9: Event yields for the evjj analysis for j6 = 0.5 for all values of Mlq. Only statistical 
uncertainties are reported, except in the "All background" column, where systematic uncer¬ 
tainties are also reported. 


MLoiGeV] 

1 Signal 

W+jets 


tt 


Multijet 

VV) Z, single t 

All background 

Data 

300 

4642 

± 

50 

822 

± 22 

1191 

± 

12 

117.9 

± 

1.5 

210.5 

± 

7.7 

2342 

± 27 

’ ± 343 

2455 

350 

2112 

± 

21 

276 

± 15 

441.4 

± 

7.2 

59.11 

± 

0.97 

102.1 

± 

5.4 

879 

± 17 

’ ± 127 

908 

400 

945.8 

± 

9.3 

110.4 

± 7.8 

184.2 

± 

4.7 

32.88 

± 

0.69 

51.5 

± 

3.8 

379.0 

± 9.' 

5 ± 53.2 

413 

450 

457.5 

± 

4.5 

53.1 

± 5.8 

74.7 

± 

3.0 

14.13 

± 

0.42 

25.7 

± 

2.7 

167.6 

± l: 

1 ± 22.2 

192 

500 

226.7 

± 

2.2 

20.5 

± 3.3 

34.4 

± 

2.0 

7.76 

± 

0.30 

15.3 

± 

2.1 

78.0 

± 4.' 

4 ± 10.1 

83 

550 

118.2 

± 

1.2 

8.6 

± 1.8 

14.9 

± 

1.4 

3.89 

± 

0.21 

7.8 

± 

1.6 

35.4 

± 2.1 

8 ± 4.5 

44 

600 

64.65 

± 

0.64 

2.3 

± 1.0 

7.08 

± 

0.93 

2.29 

± 

0.17 

4.6 

± 

1.2 

16.3 

± 1.1 

8 ± 2.1 

28 

650 

36.25 

± 

0.36 

0.41 

± 0.29 

3.82 

± 

0.70 

1.18 

± 

0.12 

2.13 

± 

0.92 

7.5 

± 1.: 

2 ± 1.1 

18 

700 

21.18 

± 

0.21 

0.41 

± 0.29 

2.61 

± 

0.60 

0.85 

± 

0.10 

0.58 

± 

0.24 

4.45 

± 0.71 ± 0.76 

6 

750 

12.56 

± 

0.12 

0.00 

-F0.94 

-n 

1.75 

± 

0.47 

0.510 

± 

0.091 

0.27 

± 

0.15 

2.54 

+1.07 

-0.50 

± 0.50 

4 

800 

7.412 

± 

0.073 

0.00 

-F0.94 

-0 

1.10 

± 

0.37 

0.317 

± 

0.067 

0.27 

± 

0.15 

1.70 

+1.02 

0.41 

± 0.31 

3 

850 

4.591 

± 

0.045 

0.00 

-F0.94 

-0 

0.90 

± 

0.34 

0.117 

± 

0.029 

0.140 

± 

0.087 

1.15 

+1.00 

-0.35 

± 0.24 

2 

900 

2.853 

± 

0.028 

0.00 

-F0.94 

-0 

0.37 

± 

0.21 

0.076 

± 

0.024 

0.084 

± 

0.069 

0.53 

+0.97 

-0.22 

± 0.10 

1 

950 

1.791 

± 

0.017 

0.00 

-F0.94 

-n 

0.37 

± 

0.21 

0.069 

± 

0.023 

0.084 

± 

0.069 

0.52 

+0.97 

-0.22 

± 0.10 

1 

1000 

1.272 

± 

0.011 

0.00 

-F0.94 

-0 

0.37 

± 

0.21 

0.069 

± 

0.023 

0.084 

± 

0.069 

0.52 

+0.97 

-0.22 

± 0.10 

1 

1050 

0.8788 

± 

0.0074 

0.00 

-F0.94 

-0 

0.37 

± 

0.21 

0.069 

± 

0.023 

0.084 

± 

0.069 

0.52 

+0.97 

-0.22 

± 0.10 

1 

1100 

0.6063 

± 

0.0049 

0.00 

-F0.94 

-n 

0.37 

± 

0.21 

0.069 

± 

0.023 

0.084 

± 

0.069 

0.52 

+0.97 

-0.22 

± 0.10 

1 

1150 

0.4196 

± 

0.0032 

0.00 

-F0.94 

-0 

0.37 

± 

0.21 

0.069 

± 

0.023 

0.084 

± 

0.069 

0.52 

+0.97 

-0.22 

± 0.10 

1 

1200 

0.2894 

± 

0.0021 

0.00 

-F0.94 

-0 

0.37 

± 

0.21 

0.069 

± 

0.023 

0.084 

± 

0.069 

0.52 

+0.97 

-0.22 

± 0.10 

1 
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7 Results 


19.7fb ' (8 TeV) 




Mei [GeV] 


Figure 9: The distribution for the eejj channel (left) and the Mq distribution for the 

evjj channel (right) after the selection criteria optimized for a LQ mass of 650 GeV have been 
applied. The dark shaded region indicates the statistical and systematic uncertainty in the total 
background prediction. The signal corresponds to a LQ mass of 650 GeV and fi = 0.075. The 
signal is multiplied by a factor of ten in the left plot. In the case of the eejj analysis, less than 
one signal event is expected to pass the selection. The horizontal lines on the data points show 
the variable bin width. 


The theoretical cross sections are represented as the central values with a band indicating the 
sum in quadrature of the PDF uncertainty and the uncertainty associated with the choice of 
factorization/renormalization scale. The latter is estimated from the observed effect of varying 
the scale between half and twice the LQ mass. 


By comparing the observed upper limit with the theoretical cross section values, first gener¬ 
ation scalar LQ with masses less than 1010 (850) GeV are excluded with the assumption that 
j6 = 1 (0.5). This is to be compared with median expected limits of 1030 (890) GeV. Similarly, 
second generation scalar LQ with masses less than 1080 (760) GeV are excluded with the same 
assumptions on /3, to be compared with median expected limits of 1050 (820) GeV. 


The combination of the ££]] and £v']i channels, shown in Fig. 12 excludes LQ masses as a func¬ 
tion of /3 using the intersection of the theoretical cross section central value and the excluded 
cross section. The combination can improve the mass exclusion reach for values of /3 < 1. 
Using the combined channels, second generation scalar LQ with masses less than 800 GeV are 
excluded for j6 = 0.5, compared with an expected limit of 910 GeV. In the case of first genera¬ 
tion LQ, the combination does not lead to a change in the observed limit for j6 = 0.5. 


The broad excess in the eejj and ei/jj channels is most significant for the final selection optimized 
for a LQ mass of 650 GeV, but has kinematic distributions that do not match those expected for 
a LQ hypothesis of that mass. Figure [T^ shows that the presence of the excess does reduce 
the exclusion power of the analysis at small values of j6 (< 0.15) for the selections optimized 
for LQ masses around 650 GeV. The exclusion limit for this region of the parameter space is 
dominated by the evjj channel. 


7.3 Additional interpretations 

Vector LQ signal samples were simulated with GalcHep at the values of LQ mass detailed 
in Section for the four scenarios of anomalous couplings described in Section The cross 
sections for pair production of vector LQs are larger than the ones for the pair production of 
scalar LQs, therefore we expect a higher reach in the Mlq exclusion limits. The cross sections 
for vector LQs have been calculated only at the LO level. We assume that the ratios of NLO to 
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Figure 10: Frame on left (right): the expected and observed upper limits at 95% CL on the LQ 
pair production cross section times (2/3(1 — /I)) as a function of the first generation LQ mass 
obtained with the eejj (evjj) analysis. The expected limits and uncertainty bands represent the 
median expected limits and the 68% and 95% confidence intervals. The left shaded regions 
are excluded by Ref. 177) and the middle shaded regions are excluded by Ref. Il24l . The right 
shaded region is excluded by the analysis presented in this paper. The CTt^eory curves and their 
bands represent, respectively, the theoretical scalar LQ pair production cross section and the 
uncertainties due to the choice of PDF and renormalization/factorization scales. 
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Figure 11: Frame on left (right): the expected and observed upper limits at 95% CL on the LQ 
pair production cross section times /3^ (2/3(1 — /3)) as a function of the second generation LQ 
mass obtained with the ijiV]]) analysis. The expected limits and uncertainty bands represent 
the median expected limits and the 68% and 95% confidence intervals. The left shaded regions 
are excluded by Ref. [78) and the middle shaded regions are excluded by Ref. Il24l . The right 
shaded region is excluded by the analysis presented in this paper. The cztheory curves and their 
bands represent, respectively, the theoretical scalar LQ pair production cross section and the 
uncertainties due to the choice of PDF and renormalization/factorization scales. 


LO cross sections for the case of vector LQs are the same as the corresponding ratios for scalar 
LQs, which vary from 1.62-4.03 over the 300-1800 GeV mass range Ill9ll . In fact, the ratios of 
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Figure 12: The expected and observed exclusion limits at 95% CL on the first (left) and second 
(right) generation scalar LQ hypothesis in the (5 versus LQ mass plane using the central value 
of signal cross section for the individual and £i/jj channels and their combination. The 
expected limits and uncertainty bands represent the median expected limits and the 68% and 
95% confidence intervals. Solid lines represent the observed limits in each channel, and dashed 
lines represent the expected limits. 


the NLO K-factors for scalar LQ pair production vs. vector LQ pair production are expected 
to be very similar to the analogous ratios for single LQ production, which have recently been 
published Il79l . Therefore, the limits we obtain by applying the scalar LQ K-factors to the 
vector LQ LO theoretical curves to obtain predictions for the NLO cross sections are expected 
to be conservative. The distributions of the kinematic variables for scalar and vector LQs are 
sufficiently similar that the same event selections and final optimization thresholds can be used 
for both analyses. It is found that the cross section limits determined using the MC scenario 
agree within uncertainties with the YM, MM, and AM coupling scenarios. Thus, it is sufficient 
to overlay the theoretical cross section curves for all vector LQ scenarios with the limit curve 
calculated using the MC scenario. 

Figurej^shows the experimental limits along with the four theoretical vector LQ cross sections 
for the eejj (evjj) channel for /3 = 1 (0.5). The experimental results yield a 95% CL upper limit 
exclusion of masses less than 1470 (1360) GeV assuming YM couplings, 1270 (1160) GeV for the 
MG couplings scenario, 1660 (1560) GeV for the MM couplings scenario, and 1150 (1050) GeV 
for the AM scenario. The increased energy and luminosity of the LHG results in considerably 
improved limits compared to the ones determined by the DO experiment at the Tevatron Il35l . 
which excluded leptoquark masses less than 340 (315) GeV for the case of YM couplings. 


Experimental limits along with the four theoretical vector LQ cross sections for the ^^jj (/^vjj) 
channel for j6 = 1 (0.5) are shown in Fig. 14 on the left (right). In the fiji]] ijiV]]) channel, the 
experimental results yield a 95% GL upper limit exclusion of masses less than 1530 (1280) GeV 
assuming YM couplings, 1330 (1070) GeV for the MG scenario, 1720 (1480) GeV for the MM 
couplings scenario, and 1200 (980) GeV for the AM couplings scenario. These are the most 
stringent limits to date on second-generation vector LQ production. 


The data have also been compared with an RPV SUSY model described in Ref. USOl . This 
model predicts light top squarks that decay to a lepton and quark through an R-parity violat- 
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Figure 13: Frame on left (right): the expected and observed upper limits at 95% CL on the 
vector leptoquark pair production cross section times (2/1(1 — j6)) as a function of the first 
generation vector leptoquark mass, obtained with the eejj (ei/jj) analysis for the four coupling 
scenarios (MC, YM, MM, and AM). The expected limits and uncertainty bands represent the 
median expected limits and the 68% and 95% confidence intervals using the MC scenario. Be¬ 
cause of the kinematic similarity between the MC scenario and the other coupling scenarios, 
cross section limits are found to be the same within the uncertainties. 
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Figure 14: Frame on left (right): the expected and observed upper limits at 95% CL on the vec¬ 
tor leptoquark pair production cross section times (2/1(1 — /!)) as a function of the second 
generation vector leptoquark mass, obtained with the /^f/jj (/^vjj) analysis for the four coupling 
scenarios (MC, YM, MM, and AM). The expected limits and uncertainty bands represent the 
median expected limits and the 68% and 95% confidence intervals using the MC scenario. Be¬ 
cause of the kinematic similarity between the MC scenario and the other coupling scenarios, 
cross section limits are found to be the same within the uncertainties. 
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8 Summary 


ing top squark-lepton-quark vertex (A') operator. The A '^32 (A 232 ) operator refers to top squark 
decay to one electron (muon) and one light-flavor quark. In the case of direct top squark de¬ 
cay, this model is kinematically similar to LQ production, and the limits already described for 
jS = 1 scalar LQs can be applied simply by scaling for the small difference in production cross 
sections between top squarks and LQs. 

It is interesting to consider the case where top squark decay is mediated by a Higgsino with a 
mass A4g = — 100 GeV with a 100% branching fraction, as shown in Fig.|^ Because of higher 

jet multiplicity and hence softer kinematic spectra, the optimization selections described in 
Section |4T] are shifted such that for a given top squark mass, the selections used correspond to 
a LQ mass lower by 100 GeV, determined by optimizing the expected limits. The experimental 
limits along with the theoretical top squark pair production cross sections for the eejj 
channel are shown in Fig. on the left (right). Assuming this model, the experimental results 
yield a 95% GL observed upper limit exclusion of top squark masses less than 710 GeV in the 
first generation A 332 model, compared with a median expected limit of 840 GeV. The second 
generation A 232 model yields an observed exclusion of top squark masses less than 860 GeV, 
compared with a median expected limit of 880 GeV. These are the first experimental limits to 
date on A 332 and A 232 RPV SUSY top squark decays. 
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Figure 15: Frame on left (right): the expected and observed upper limits at 95% GL on the top 
squark pair-production cross section for a Higgsino-mediated RPV SUSY model in the eejj (f/;^jj) 
+ 4 b quark final state as a function of the top squark mass, obtained with the eejj analysis. 
The expected limits and uncertainty bands represent the median expected limits and the 68 % 
and 95% confidence intervals. The chthgory curves and their bands represent, respectively, the 
theoretical top squark pair production cross section and the uncertainties due to the choice of 
PDF and renormalization/factorization scales. 


8 Summary 

A search has been conducted for pair production of first- and second-generation scalar lep- 
toquarks in final states with either two electrons (or two muons) and two jets, or with one 
electron (or muon), significant missing transverse energy, and two jets, using 8 TeV proton- 
proton collision data corresponding to an integrated luminosity of 19.7 fb~^. The results are 
also interpreted in the context of models of vector leptoquark pair production and of R-parity 
violating supersymmetric models with similar final state signatures. 
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The selection criteria used for all the searches are optimized for each scalar leptoquark signal 
mass hypothesis. In the first generation eejj (ei/jj) channel, a broad 2.3 (2.6) standard deviation 
excess is observed in the final selection optimized for leptoquarks with a mass of 650 GeV. The 
excess does not peak in the Mq distributions, as a leptoquark signal would, but does weaken 
the upper limit that can be set on the production cross section for leptoquark masses of about 
650 GeV and values of /I < 0.15. Limits are placed with 95% GL on first-generation scalar lep¬ 
toquarks with masses less than 1010 (850) GeV, assuming j 6 = 1.0 (0.5). This is to be compared 
with the expected 95% GL exclusions of 1030 (890) GeV. In the secondgeneration leptoquark 
search the number of observed candidates for each mass hypothesis agrees within uncertain¬ 
ties with the number of expected standard model background events. Second-generation scalar 
leptoquarks are excluded at 95% GL with masses below 1080 (800) GeV for /3 = 1.0 (0.5). This 
is to be compared with a median expected limit of 1050 (910) GeV. These results for pair pro¬ 
duction of scalar leptoquarks are closely comparable to those of Ref. II23II . 

Limits are set on four coupling scenarios for vector leptoquarks, and for the eejj (evjj) channel 
yield 95% GL upper limit exclusions of masses in the range of 1150 — 1660 (1050 — 1560) GeV. In 
the }inii (f/vjj) channel, the experimental results yield 95% GL upper limit exclusions of masses 
in the range of 1200 — 1720 (980 — 1480) GeV. These represent the most stringent limits on 
vector LQ production to date. 

Limits are also set for top squark production in an R-parity violating supersymmetric model 
via the Aj 32 or A 232 operators. For direct top squark decay, the scalar LQ limits can be applied 
directly. Interpretation is also made in Higgsino-mediated top squark decay, where the exper¬ 
imental results yield a 95% GL observed upper limit exclusion of top squark masses less than 
710 GeV in the first generation A 332 model, compared with a median expected limit of 840 GeV. 
The second generation A 232 model yields an observed exclusion of top squark masses less than 
860 GeV, compared with a median expected limit of 880 GeV. These represent the most strin¬ 
gent experimental limits to date on A 332 and A 232 RPV SUSY top squark decays and the first 
experimental limits on the Higgsino-mediated decays. 
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